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Sources of internal gravity waves (IGW) in the upper atmosphere are 
assumed to be meteorological processes in the troposhere. These sources 
are vertically and horizontally inhomogeneous and time-dependent. In order 
to describe the IGW propagation from such sources, a numerical solution of 
a system of hydrodynamical equations is required. In addition, it is 
necessary to take into account the influence of the altitude-latitude 
inhomogeneity of the temperature and wind fields on the IGW propagation as 
well as the processes of dissipation. 

This paper proposes an algorithm for numerical modelling of the IGW 
propagation over a limited area from tropospheric local sources to the 
upper atmosphere. The algorithm takes into account all the above mentioned 
features. A spectral-grid method is used with the expansion of wave fields 
into the Fourier series over longitude. The upper limit conditions have 
been obtained from the requirement of a limited energy dissipation rate in 
an atmospheric column. The no slip (zero velocity) condition has been used 
at the Earth's surface. 

Therefore, for the analysis of the vertical structure of steady-state 
oscillations, the IGW source is given in the equation for divergence as an 
inhomogeneous term 


0 = f • exp [ - 
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where f characterizes the source intensity, d its vertical extension, z is 
the height of the maximum, k and a are the horizontal wave number and the 
IGW frequency, respectively. The values of Z q = 10 km and d j=^5 _lpn have 
been used in the calculations, as well as the value of f = 10 s . The 
results obtained (see below) show that even this moderate value of f 
provides correct orders of magnitude of the IGW amplitudes in the upper 
atmosphere . 


Wave equations are integrated over time after switching on the source 
(1) at a movement t = 0 prior to the transition of the solution to the 
stationary regime. For complete stabilization of the wave oscillations 
regime, a time At-(60-80)T is required (t is the IGW period). 


Fig. 1 shows profiles of the variations of wave energy dissipation e ; 
of the transmission of wave energy to the background wind e and the wave 
acceleration of the background flux a as well as the rate of dissipation of 


wave energy into heat e 
lower thermosphere. 


All curves have maxima at the heights of the 
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Since for the stabilization of oscillations in the atmosphere time 
intervals of tens and hundreds of hours are required, as was mentioned 
above, real meteorological sources in the atmosphere cannot, usually, 

remain unchanged during such long time intervals. For modelling such a 

source the above-mentioned inhomogeneous term in that equation for 
divergence has been taken as 

<f>(x.z.t) = f exp [ - ( t - t ) 2 /t 2 - (z - z ) 2 /d~ - ik x], (2) 

1 1 O O O X 

where t is the scale of variation of the wave source in time and t is the 
time of°its maximum. It can be considered that the source (2) generates a 
packet of I GW with the frequencies -2 tt/t <u>< 2tt/t o which, with a given k^, 

are propagated in both directions of the axis with horizontal phase 

velocities |C | < C = 2 tt/(t k ). Therefore, the values t and C have 

been used as parameters characterizing non-stationary sources’. In ("2) the 
following parameters for a tropospheric source are used: 

f = 10^s 2 , z = 10 km, d = 5 km 
o 

The results of calculations of the time variation of the integral over 
an atmospheric column of the wave energy E tota i show that after the 
"switching on" of the source the wave energy quickly reaches its maximum 
values and then a slow process of attenuation begins which continues for 
(20-40) t . This can be accounted for by the fact that the main part of 
wave energy is concentrated in the tropo-strato-mesophere and by the 
relatively low speed of its propagation to the thermosphere where 
dissipation of waves mostly takes place. The time lag of the wave energy 
maxima in relation to the t /5 moment t for E . is from t / 5 to t and 

for local values of t E °increases for up to a C"4-8) t at 100-150 km 
, . , o o 

heights . 

Fig. 2 shows the vertical profiles of the IGW parameters for January 
at various times for sources with t =6 hrs and different C . Analysis of 
Fig. 2 shows that after a time period (4-8) t following the source 
function maximum a comparatively stable profile of E is formed with two 
maxima: in the troposphere and at the height of about 75 km. These maxima 

coincide with the regions of minima gH where g is the static stability 
parameter. The formation of the mesopheric waveguide is additionally 
influenced by the background wind profile which creates an increase of the 
wave flux F at the decrease of A near the 75 km height. The formation of 
the waveguide leads to the lowering of the maxima of wave amplitudes and 
e , e, and a in Fig. 2 in comparison with the stabilized regime. Curves 
K* and A in Fig. 2 show that in the initial period the main reflection of 
wave energy takes place at 80-90 km due to the effect of the temperature 
profile and the maximum of turbulent dissipation in the lower thermosphere. 
As time passes, this barrier grows more transparent and another area of 
increased reflection is formed near the jet stream at a height of about 50 
km. 


In most cases the background wind in the strato-mesophere leads to the 
enhancement of the harmonics of IGW coming from the opposite direction in 
comparison with those propagating in the same direction, and the resulting 
propagation of the IGW packets is directed towards the background flux. 
This determines the signs of accelerations a of the background flux by 
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Vertical profiles of IGW parameters for January at various times 

Dashed, solid and dashed-dotted lines correspond to c = 90, 150 

~m 

and 250 mps. 


ORIGINAL PAGE IS 
OF POOR QUALITY 






171 


waves. In winter a<o (directed westwards); in summer, at its beginning and 
at heights less than 100 km, a>o (directed eastwards) although after a time 
(8-12) t after t and above 100 km a<o occur. According to the 
present-day understanding (HOLTON, 1983), the effect of IGW drag is an 
essential mechanism of formation of the middle atmosphere circulation. The 
calculations presented in this paper, according to which the filter formed 
by the background wind, for the most part, transmits IGW propagating from 
the opposite direction, account for the mechanism of appearance of dragging 
wave accelerations. 

Fig. 3 pre^ei^ts the dependence on = o/k^ of the normalized wave 
energy E . C /C integrated over an atmospheric column, where C = 150 
mps, and on the densities of wave energy E at different height levels. The 
solid curves in Fig. 3 reveal, for a windless atmosphere, the resonance 
peaks at the values C = 310, 255 and 180 mps which are very close to the 
values of atmospheric oscillations discussed by DIKIY (1969). The 
influence of the background wind leads to a bias in the resonance maxima 
for IGW coming from the direction opposite to the wind to C = 295, 215 and 
125 mps, and for IGW propagating in the direction of the wind to C = 335, 
300 and 235 mps, respectively. The bias of the resonance maxima in Fig. 3 
is accompanied by a decrease in their values for IGW coming from the 
opposite direction and by an increase for those propagating in the same 
direction. Analysis of Fig. 3 (b through f) shows that, in the lower 
layers, the first resonance energy maximum prevails which corresponds to 
C ~300-315 mps. However, as the height increases it is attenuated and at 
z~100-150 km the second resonance maximum prevails. 
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Fig. 3 Dependence on c = o/k of E (frame a) and E at the height 
levels 1, 10, 58, 100 and 150 km, (b) through (f). 



